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Abstract. We consider the mean first passage time of a random walker moving in 

a potential landscape on a finite interval, starting and end points being at different 
C^ \ potentials. From analytical calculations and Monte Carlo simulations we demonstrate 

that the mean first passage time for a piecewise linear curve between these two points 
' , is minimised by introduction of a potential barrier. Due to thermal fluctuations this 

^H I barrier may be crossed. It turns out that the corresponding expense for this activation 

O ' is less severe than the gain from an increased slope towards the end point. In particular, 

the resulting mean first passage time is shorter than for a linear potential drop between 

the two points. 
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1. Introduction 

In classical mechanics, Bernoulli's 1696 brachistochrone problem addresses the curve 
between two points that is covered by a point particle in the least time, under the 
influence of gravity. If the particle starts at rest the brachistochrone curve is a cycloid. 
Steeper at first, the particle is accelerated, keeping its momentum in absence of friction. 
In particular at no point along this curve the particle elevation is higher than that of the 
starting point, for reasons of energy conservation. An overdamped, diffusing particle 
may appear to behave classically: driven by a constant external force the mean first 
passage time (MFPT) T from one point to another along the direction of the force 
equals L/V, the ratio of distance L versus the particle velocity V [1]. However, as the 
diffusing particle is coupled to a heat bath, thermal fluctuations may lift it across a 
potential barrier. At the same time, the overdamping does not allow the particle to 
take along its momentum. To minimise the MFPT one would thus naively expect that 
the particle should constantly move downhill. As we are going to show here for the case 
of a piecewise linear potential, it is indeed beneficial for the MFPT if the particle first 
crosses a potential barrier, that is, the particle initially moves uphill. As a consequence 
the following downhill slope becomes steeper, leading to a smaller overall MFPT. 

Generally, the question of the interplay between potential landscape and diffusion 
properties is of great interest, resulting in often surprising behaviour such as giant 
diffusivity [2]. But which shape of the potential should one choose in order to optimise 
the escape time on an interval? A large number of previous studies were concerned with 
problems of the escape from a potential well [3], following Kramers' classical work |3]. 
Optimisation of the escape time may involve phenomena such as resonant activation 
[5] . One of the simplest models for a potential landscape is a piecewise linear potential 
(Fig. [1]). Only recently it was realised that an asymmetry in this kind of potential 
is important for escape properties in resonant activation [HI [7]. The asymmetry of 
the potential also plays a crucial role in systems with periodic potentials relevant to 
molecular motor models |8l|9l[10], or for molecular shuttles in suprachemical compounds 
|11] . However, to the best of our knowledge the role of asymmetry on the MFPT for a 
static potential as displayed in Fig. [1] has not been discussed. 

2. Mean first passage time minimisation 

We consider a particle diffusing from the starting point O at x = 0, to point X located at 
xx = 1, in the piecewise linear potential going through point A at xa- This situation is 
sketched in Fig. [1] The values of the potential in these points are Uq, Ua, and Ux = 0, 
without loss of generality. At the starting point O we impose a reflecting boundary 
condition while at the end point X we apply an absorbing boundary condition for the 
calculation of the MFPT. The question we pursue is: which shape of the piecewise linear 
potential minimises the MFPT from O to X7 

The MFPT for the piecewise linear potential with bias velocities vi (on < a; < xa) 
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Figure 1. Scheme of the piecewise hnear potential (blue line) considered here. The 
particle is initially placed at point O (at x — 0), in which we impose a reflecting 
boundary condition. The end point is X, and we choose xx = 1- At the turnover 
point A the slope of the potential changes. v\ and V2 are the drift velocities on the 
two linear slopes (v\ < and t;2 > 0). The red line shows the inversely symmetrical 
potential resulting in the same MFPT (see text). 



and V2 (on xa < x < 1) on the unit interval, shown in Fig.lH readily obtains analytically 
[1|,[T2]. A unit current j(0,t) = 6(t) is injected at x = 0, and the output is calculated 
from the solution of the Fokker-Plank equation, 
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where U'{x) is the derivative of the external potential. Moreover m is the particle 
mass, r} the friction experienced by the particle, and D is its diffusion constant. For 
the gravitational potential U{x) = mgh{x) for a particle at elevation h{x) at position 
X and with the gravitational constant g, the drift term in the Fokker-Planck equation 
becomes d/dx {gh\x)/r}) P{x, t). The ratio g/rj has the dimension of a velocity, so that 
the Fokker-Planck equation may be rewritten in the form 
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with piecewise constant drift velocity Vi, where i = 1,2 denotes the two domains with 
piecewise linear potential. Note the sign of the drift velocity: an increase of the potential 
causes a drift to the left, and vice versa. The reflecting and absorbing boundary 
conditions at x = and x = 1, respectively, read ^\ and P(l,t) = 0. Requiring 
continuity of the distribution P and the probability flux at point A, the MFPT yields 
in the form [1] 

D 



T 



VxV-i 



(1 



-vvxaI'D 



)(1 



-~V2{\-XA)|D^ 



Xa 1 

+ — + - 



Xa D (e' 



v\xa/D 



1) 



V2 



How a finite potential barrier decreases the mean first passage time 4 

D (q-V2{1-xa)/D _ ^\ 

+ — -2 -, (3) 

as function of xa, vi, and V2. We note that all variables occurring in Eqs. (II]) to ([3]) 
are dimensional. In what follows we measure lengths in units of cm and time in sec. 
Thus when writing L = 1 for the distance between starting and end points, this actually 
means 1cm. 

Let us study the MFPT ([3]) in detail. We first note that expression (|3]) is symmetric 
under simultaneous exchange of vi ^ t>2 and x^ ^ 1 — xa, i.e., inversion through the 
midpoint of the line connecting O and X. This inverse case corresponds to the red line 
in Fig. [H Secondly, we observe that by increasing the elevation of point A with respect 
to O and X and shifting the turnover point A towards the starting point O such that 
I^iXaI ^ 1, \v2\{'^ — Xa) ^ 1, and xa <^ 1, the MFPT ^ reduces to 

T ^ :^el"^"-^l/^ + -. (4) 

This is the sum of the MFPTs on the two subintervals. Indeed, the first term corresponds 
to the Kramers rate for crossing of a high potential barrier, see below, while the second 
term represents the MFPT at constant drift V2 over the unit distance. Result (jl]) 
demonstrates that the the overall MFPT T as well as both individual terms are reduced 
by increase of A's elevation while keeping the product viXa constant. This is one of the 
central results of our study: the introduction of a high but narrow barrier reduces the 
MFPT. 

For the thermally activated crossing of a sufficiently high potential barrier the 
corresponding barrier crossing time was obtained by Kramers [H [13] . 

Here Xmin and Xmax denote the positions of the potential minimum (where the particle 
is initially placed) and the saddle of the potential. According to expression (^ this 
characteristic time depends both on the potential difference AU = ^/(xmax) — U^x^i^) 
and the curvature of the potential at these two points. If we imagine that we smoothen 
the piecewise linear potential around the minimum and maximum points, it becomes 
clear that for fixed AU a decrease of the distance between Xmin and Xmax implies an 
increase of the respective curvatures and thus a decrease of the barrier crossing time. 
This observation underlines that our above result for the MFPT in the piecewise linear 
potential is consistent with the physics of barrier crossing. 

What happens in the case opposite to Eq. (II]), when the two drift velocities are 
small, |fi|, \v2\ ^ 1? Expansion of Eq. (^ up to first order yields 
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Figure 2. Minimal MFPT in the piecewise linear potential for vanishing potential 
difference between starting and end point, as function of the turnover point position 
XA (blue curve) . The corresponding optimum value for the value of the potential at the 
turnover point is shown as the red line. The dashed line emanating from the turnover 
point in the schematic of the potential profile (bottom of graph) intersects the two 
curves at the associated values of MFPT and Ua- 



Here, the first term represents the MFPT of free diffusion on the unit interval. The next 
two terms are the first order corrections in Vi and V2- Depending on the actual values 
of f 1 and f 2 these terms may either lead to a decrease or increase of the MFPT. 

While the MFPT can be arbitrarily reduced by increasing vi (and thus also V2) and 
simultaneously decreasing the position xa of the turnover point, a finite potential barrier 
may still reduce the MFPT. We analyse the three possible, different cases in Figs. [JMl 
Starting with the case when starting and end points are at the same potential level, in 
Fig. [2] we show the minimal value for the MFPT (jS]) together with the corresponding 
optimal value for the potential at A, Ua, as function of the position xa of the turnover 
point. This minimisation was performed numerically with Mathematica. We see that 
the largest value of the MFPT is obtained when the turnover point is located in the 
middle of the interval at xa = 0.5. In this special case the optimum is reached in absence 
of a potential barrier {Ua = 0), i.e., for unbiased diffusion. Away from the midpoint, 
the MFPT appears dramatically reduced. For xa -^ and a;^ — )■ 1, the fastest MFPT 
is obtained when the potential diverges, Ua — ^ ±00. Notice the symmetries of both 
MFPT and profile of optimal turnover points with respect to the midpoint, xa = 0.5. 

For the case of very asymmetric positions of turnover points xa — > 0, the 
optimal value for the drift Vi can be computed analytically, if the potential difference 
Af/ = Vix + f 2(1 — x) and xa are fixed. Expansion of expression ([3]) as a series for small 
Xa leads to the first order approximation 
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Figure 3. Minimal MFPT and corresponding height of the potential at the turnover 
point A as function of the position xa, in the case when the potential difference between 
starting and end points is 10 fcsT. 



where Q± = exp{—AU/D) ± 1. Here the first two terms are the MFPT for a uniform 
hnear bias with potential difference AU. The third term is the correction hnear in xa- 
Analysing its form shows that an increase of the height of the turnover point (i.e., an 
increase of |fi|) always leads to a decrease of the MFPT if AU is positive. For the 
optimal slope vi we obtain the approximate expression 



Vl 



AU {^n+ + 2n^) (2 - AxA - ^xa) 
'2xA (Qn_ + f?e-^ + 2^(1 + 2e-^) 



(8) 



In the range of small xa and AU > all terms in the brackets are positive. Hence, 
expression ([H]) proves analytically that in this case a barrier indeed optimises the MFPT. 
Note that the numerical accuracy of this approximation is actually not too good. In 
order to reproduce the functional behaviour over a longer range of xa higher order terms 
need to be considered. 

For the case when the starting point is higher than the end point, the result for 
the minimal MFPT is displayed in Fig. [3l Here the MFPT shows an extended plateau 
around xa = 0.5. Exactly at this midpoint the minimum MFPT corresponds to the 
naively expected case of a constant slope from starting to end point. For xa closer to 
zero the MFPT again drops down to zero while the value of the potential at the turnover 
point diverges. Both curves for the MFPT and the potential at the turnover point are 
again symmetric with respect to the midpoint. In contrast to Fig. O however, the curve 
for the MFPT is not symmetric around the zero-line of the potential. 

For completeness we consider the case when the end point is elevated with respect 
to the starting point. While a classical particle would never reach this end point, 
a thermally driven particle may gain the necessary energy from the heat bath. The 
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Figure 4. Minimal MFPT and associated turnover potential for the case when the 
potential difference between starting and finishing points is — lO/csT (end point is 
higher than the starting point). 



corresponding optimal potential of the turnover point in the piecewise linear potential 
and the associated MFPT are shown in Fig. HI It turns out to be beneficial when an 
initial barrier exists whose height exceeds the overall potential difference |A[/| between 
starting and end point, such that the drift velocity v^ is positive. 

Let us compare the minimal MFPT in the three cases of positive, zero, and 
negative potential difference between the initial and end points of our setup, for 
XAi = 10-3 and xa2 = 0.5 (i.e., the longest MFPT). For AU = lOkeT (Fig. [3]) the 
ratio T{xai) : T{xa2) ~ 0.7, for AU = (Fig. [2]) it is T{xai) : T{xa2) ~ 0.29, and for 
AU = -lOkeT (Fig. Hj) we find T{xai) : T{xa2) ~ 0.01. Thus, the introduction of a 
potential barrier or kink has indeed the largest effect on the MFPT when the end point 
has a higher energy. That is, when it is harder to reach the end point energetically, the 
benefit from a potential turnover is larger. This is the second central result of our study. 

We simulated the Brownian motion of a particle in a piecewise linear potential with 
a Monte Carlo approach, based on the Metropolis algorithm: If the potential difference 
6U between current and potential new position is positive, 6U > 0, then the step is 
accepted with probability exp {—6U/[kBTM]), where ksTM is a measure of temperature. 
Otherwise the step is immediately accepted. 

Comparison with the analytical results was achieved by consideration of the 
continuum limit of a discrete biased random walk on a lattice. The probability 
distribution of jumps of length £, p(£), is defined by the Fokker-Planck equation [14j 

(9c(x, t) 
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where it is assumed that the lattice spacing and time step are infinitely small: A — )■ 0, 
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Figure 5. Comparison of Monte Carlo simulations (squares) with the analytical result 
from Eq. ([3]) shown as the full line. The lattice size is N ^ 1,001, and the number of 
runs is 100,000. 
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In the case we considered, the values of diffusion constants and the slopes in continuum 
limit are 
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2ArV 2xAN^kBTMr' 2{1 - XA)N^kBTMr' 

where N is the lattice size and xa the position of the turnover point. 

The simulations demonstrate excellent agreement with our analytical results. We 
show the comparison between simulations and Eq. ([3]) for the case AU = for xa = 0.1, 
kj^TM = 1, and A^ = 1001 in Fig. E 



3. Discussion 

On a flat potential landscape significant progress has been achieved in the theory of 
MFPTs on arbitrary, finite domains [15]. In particular, the role of compact versus non- 
compact explorations has been revealed in generality [TB]. Much less is known about 
MFPT properties in potential landscapes. 

We analysed the value of the MFPT in a finite interval for a piecewise linear 
potential, finding that the introduction of a barrier reduces the MFPT. In the ideal case 
when the barrier height is unlimited the MFPT can be reduced arbitrarily. This a priori 
surprising results were shown to be in line with physical arguments such as Kramers 
escape theory, and may be of interest in the design of potential energy landscapes, 
for instance, for functional molecules (molecular shuttles), or for molecular motors. 
Conversely, our results may shed new light on the role of barriers in known landscapes, 



How a finite potential barrier decreases the mean first passage time 



v,= j_ 




Figure 6. The MFPT in case (a) is larger than in case (b) as long as Ua is fixed. 
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Figure 7. MFPT for the case of fixed potential Ua corresponding to Fig. IH^b) as 
function of the position xa of the turnover point. Here, AU = 5k bT. 



for instance, in the folding landscape of proteins. Indeed, it was shown in Ref. [IT] that 
intermediate barriers of height > IfcsT increase the folding rate of proteins |T7]. 

What happens if the height of the potential barrier is limited? Consider the 
situation sketched in Fig. [61 If the values of U at the starting and end points of the 
interval are fixed and the height of the potential in point A is fixed, it is clear that in 
case (a) the MFPT is higher than in case (b). This changes considerably the answer to 
the MFPT minimisation task. Starting with a horizontal slope we could still imagine 
that a shift of the turnover point A may optimise the MFPT: if shifted to the right we 
have an increase for the time to reach A but a gain from an increased drift velocity V2- 
Variation of xa in this case leads to the dependence shown in Fig. [71 At the right end of 
the interval between starting and end points the behaviour tends to the value T = 0.5, 
corresponding to unbiased diffusion. The gain at the optimum value for xa in this case 
is in fact only a few per cent, compared to the case of a linear potential drop {xa = 0). 

All results presented above demonstrate the critical importance of the asymmetry 
of the potential barrier for optimisation of the MFPT. This gain rests on the significant 
facilitation of the passage on the long easy slope which overcompensates losses for 
crossing of the barrier. The result ([3|) allows the adjustment of the MFPT to any 
finite value, including infinitely large and infinitely small times. However, if one wants 
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to decrease the MFPT to some specific, small value, this result also shows that, to 
compensate an increase in barrier height, a substantial reduction of the position xa of 
the turnover point is required. 

In classical mechanics the cycloid is the optimal curve for a point particle in absence 
of friction: after an initial steep descent, i.e., high acceleration, the momentum of the 
particle carries on. For a diffusing, overdamped particle in the case of piecewise linear 
potential provide the answer is qualitatively the opposite: in order to minimise the 
MFPT there should be a steep and short ascent. 

It will be interesting to consider more complex shapes of the potential, in particular, 
the case of multiple barriers as mentioned in the context of protein folding [T7] . 
Moreover, numerical analysis of the first passage distribution associated with the process 
considered herein will be of interest, as well as the consideration of the full motion 
including inertia effects. Although it is possible to optimise the potential by trial and 
error for a fixed set of potential shapes, the question about whether the optimisation 
algorithm exists in generality, remains to be investigated. Another interesting question 
is whether similar results could be obtained under anomalous diffusion conditions [18]. 
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